Phase separation on a hyperbolic lattice 
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We report a preliminary numerical study by kinetic Monte Carlo simulation of the dynamics 
of phase separation following a quench from high to low temperature in a system with a single, 
conserved, scalar order parameter (a kinetic Ising ferromagnet) confined to a hyperbolic lattice. 
The results are compared with simulations of the same system on two different, Euclidean lattices, 
in which cases we observe power-law domain growth with an exponent near the theoretically known 
value of 1/3. For the hyperbolic lattice we observe much slower domain growth, consistent to within 
our current accuracy with power-law growth with a much smaller exponent near 0.13. The paper 
also includes a brief introduction to non-Euclidean lattices and their mapping to the Euclidean 
CN plane. 
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I. INTRODUCTION 



i-G The geometry most familiar to condensed-matter physicists is the Euclidean one with its vanishing Gaussian cur- 

«"n vature [1]. The circumference and area of a circle of radius p in the Euclidean plane are the well-known power laws, 

C{p) = 2irp and A(p) ~ J^C(r)dr ~ np 2 , respectively. Almost equally familiar are elliptic or spherical surfaces, 
exemplified in the macroscopic world by the Earth's surface and in the nanoscopic world by carbon buckyballs. These 
closed surfaces have positive Gaussian curvature, k > 0. Hereafter using dimensionless units such that |k| = 1, the 
circular circumference and area are analogously given by C(p) = 2i:svo.p and A(p) = 27r(l — cosp). 
^ More exotic to most is probably the hyperbolic geometry with its negative Gaussian curvature, k < 0. In this 

case the dimensionless circular circumference and area are exponentially divergent: C(p) = 27rsinhp and A(p) — 
27r(coshp — 1). The best known example is the Minkowski metric of relativistic spacetime. However, hyperbolic 
surfaces have recently been studied in nanoscience as well, including junctions of several carbon nanotubes [5] and 

i-£j anisotropic lipid membranes [3J. Percolation on hyperbolic lattices has also been studied [I]. 

In this paper we present preliminary results on a comparison of the dynamics of pattern formation during phase 
separation (spinodal decomposition) in media confined to Euclidean and hyperbolic surfaces. As our example we use 
an S = 1/2 ferromagnetic Ising model on a regular lattice embedded in the surface, and we study the time evolution 
of the characteristic pattern length following a quench from infinite temperature to one well below the model's critical 

i— ( temperature. 

The rest of the paper is organized as follows. In Sec. [H] we describe the Poincare disk mapping used to map 
patterns on a hyperbolic surface onto a Euclidean plane. Next, in Sec. |IH| we describe the lattices generated by 
regular tesselations of Euclidean, spherical, and hyperbolic surfaces. The phenomenology of phase separation is 
briefly reviewed in Sec. IV and the methods of simulation and data analysis are discussed in Sec. [V] Our numerical 



results are presented in Sec. |VI[ and a final discussion and suggestions for future work are given in Sec. |VII 
II. POINCARE DISK MAPPING OF NON-EUCLIDEAN SURFACES 



Patterns on a non-Euclidean surface cannot be reproduced on a Euclidean surface without distortion; angular, 
linear, or both. No single rendering is ideal in every respect, as evidenced by the many different geographic map 
projections developed over the centuries. The projection of the hyperbolic plane that we use in this paper is known 
as the Poincare disk mapping and is illustrated in Fig. [11 It consists of a Euclidean plane, a unit sphere (n = +1), 
and a unit hyperboloid of revolution (n = —1) 5]. The plane forms the equatorial plane of the sphere, and the 
hyperboloid rests with its apex on the North Pole of the sphere, (0,1). A straight line through the South Pole, (0,-1), 
connects the points H on the hyperboloid and S on the sphere with their mapping D on the plane. The region of 
the equatorial plane inside the sphere is the Poincare disk. It is easily seen from Fig. [TJthat an arbitrary point H on 
the hyperboloid is mapped inside the Poincare disk, leading in the map to an exponential contraction of lengths far 
away from the apex of the hyperbola. The mapping from H to D is conformal (angle-preserving) , and geodesies on 
the hyperboloid are mapped onto circles on the disk that meet its edge at straight angles. A little calculus shows that 
the "hyperbolic radius" of H, pn, is the arc length of the geodesic from (0,1) to H, calculated with the Minkowski 
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FIG. 1: The Poincare disk mappings from a unit hyperboloid of revolution (k — —1) and a unit sphere (k — +1) to a Euclidean 
plane. The y coordinate points into the page and has been suppressed in our notation, describing the position of a point as 
(x,u). See discussion in Sec. |ll| 

metric, ds — \Jdx 2 — du 2 . (Analogously, the "spherical radius" of S, ps, is the arc length of the geodesic from (0,1) to 
S, calculated with the Euclidean metric, ds = y/dx 2 + du 2 .) Other equations relevant to the mappings are included 
inFig.[T} 

III. EUCLIDEAN, ELLIPTIC, AND HYPERBOLIC LATTICES 



To construct lattices embedded in Euclidean and non-Euclidean surfaces, we consider regular tessellations of such 
surfaces by regular polygons. A lattice created by this procedure, such that q regular p-gons meet at every lattice 
site is characterized by its Schlaffli symbol, {p, q} [I] . A lattice of finite size is often denoted by the amended Schlaffli 
symbol, {p, q, R}, where R is the number of concentric layers of p-gons surrounding the central site. The only three 
regular tessellations of the Euclidean plane are shown in Fig. [2] Any regular p-gon can be decomposed into p isosceles 
triangles that meet at its centroid. This is not only true for Euclidean lattices. An illustration for the hyperbolic case 
is shown in Fig. 3] Each triangle has apex angle 9 and basal angles 0/2. 

For the Euclidean plane, the interior angles of a triangle must always sum to it. Thus, 



2?r/p + 2n/q = n <& (p - 2)(q - 2) = 4 . 



(1) 



This proves that the only combinations of integer p and q compatible with Euclidean geometry are {3, 6}, {4, 4}, and 
{6,3}. 

Similarly, for the elliptic plane, 



2?r/p + 2n/q > it <& (p - 2)(q - 2) < 4 . 



(2) 



Again it is clear that the number of possible regular tessellations is finite. In fact, they correspond to the five Platonic 
solids, {3, 3}, {3, 4}, {3, 5}, {4, 3}, and {5, 3} |B]. 
For the hyperbolic plane, on the other hand, 



2tt/p + 2n/q < it <3> (p - 2){q - 2) > 4 . 



(3) 
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{4,4} 



{6,3} 



FIG. 2: The three regular Euclidean lattices {p, q} of q p-gons meeting at each lattice site. Here, each is shown with r — 2 
layers of polygons: {3, 6, 2}, {4, 4, 2}, and {6, 3, 2}. The {3, 6} and {6, 3} lattices are each other's duals, while {4, 4} is self-dual. 
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FIG. 3: Splitting a hyperbolic pentagon into five isosceles triangles. 



Consequently, the number of possible regular tessellations is infinite. Some examples are shown in Figs. [4] and [5] 
As the size of the lattice increases, embedding without overlaps into a three-dimensional Euclidean space becomes 
impossible. Fascinating images of models of hyperbolic planes created by crocheting can be found in Ref. [7]. 



IV. PHASE SEPARATION 



Phase separation occurs when a binary mixture is quenched from a high temperature into the phase-coexistence 
region below its critical temperature. As coherent regions of the two coexisting phases form and grow after the quench, 
the length scale characterizing the typical domain size increases algebraically with time as 



f 



(4) 



The growth exponent n depends on the symmetries and conservation laws governing the dynamics. For the case of 
two-phase coexistence and a constant ratio of the volumes of the two phases, n = 1/3. This situation is known in the 
terminology of critical dynamics as Lifshitz-Slyozov dynamics [8 j or Model B \§\. However, these results implicitly 
assume a Euclidean geometry, and we are not aware that the prediction has yet been tested in the hyperbolic case. 
A numerical test is the purpose of the work presented here. 




{7,3} {4,5} {5,4} 

FIG. 4: Some examples of the infinite number of hyperbolic lattices, projected onto the Poincare disk. The exponential length 
contraction in the projected image of large pu is clearly evident. As for the Euclidean and elliptic geometries, {p, q} and {q,p} 
are duals. 




FIG. 5: Perspective image of a small {3,7,4} lattice embedded in a three-dimensional Euclidean space. The exponential 
divergence of the hyperbolic circumference with radius prohibits such embedding without overlaps for large lattices. 



V. SIMULATION AND DATA ANALYSIS 



We consider the phase separation that occurs when a S = 1/2 Ising fcrromagnet with spins Sj = ±1 placed at 
the vertices of the {p, q} lattice is quenched from a high temperature to one well below its critical temperature. The 
Hamiltonian is given by 



H = -J'£ 



IjOj . 



(5) 



(Uj) 



Here, J > is the ferromagnetic interaction constant, and the sum runs over all nearest-neighbor pairs. The co- 
ordination number in a {p, q} lattice is q. We will use units such that Boltzmann's constant and J both equal 
unity. 

While the phase transition at the critical temperature T c for this model on Euclidean lattices belongs to the two- 
dimensional Ising universality class, on hyperbolic lattices it belongs to the mean- field universality class [TDJ [TT]. In 
either case, the value of T c increases with q. (T c = 2/ln(l + y/2) « 2.269 for {4,4}, 4/ In 3 « 3.641 for {3,6}, and 
~ 5.5 for {3, 7} |11|.) To minimize surface effects, we use periodic boundary conditions for the two Euclidean lattices. 




FIG. 6: Snapshots of the {3, 6, 60} (left) and {3, 7, 8} lattices (right) in their initial, disordered states and at times near the 
end of the simulation runs. T = 2.0, in both cases well below T c . Small thermal fluctuations in the bulk phases are noticeable 
for the {3, 6} case. 

Unfortunately we are not aware of a method for doing so in the hyperbolic case, and consequently we simulate the 
{3, 7} lattice with free boundary conditions. 

The initial state is a random distribution of up and down spins, subject only to the constraint of a vanishing order 
parameter, £3 . s i = 0. The time evolution is obtained from a kinetic Monte Carlo (MC) simulation by the order- 
parameter conserving Kawasaki dynamics |12] . This algorithm consists in randomly choosing a nearest-neighbor spin 
pair and checking if the two spins are different. If they are equal, a different pair is chosen. If the spins are different, 
they are exchanged with the Metropolis probability, 



Ji exv^i 



min[l,exp(-A£/T)] 



(6) 



where AE is the energy change that would result from a successful spin exchange. In a system consisting of N spins, 
N random choices of a spin pair constitute the MC time unit, one MC step per spin (one MCSS). Snapshots of 
the {3,6} and {3,7} lattices in their initial, disordered states and at t = 10 6 MCSS, when macroscopic bulk phase 
domains are well developed, are shown in Fig. [6] 
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FIG. 7: The spin correlation function G(r) for a {3, 7, 7} lattice at different times. The correlation length f (t) is estimated as 
the distance corresponding to the first zero crossing of G(r), as discussed in the text. 



The power-law result for the characteristic length scale given in Eq. Q assumes an isotropic system with sharp 
interfaces and no thermal fluctuations in the bulk phase regions. Neither assumption is well satisfied for discrete Ising 
models at nonzero temperature. Care must therefore be exercised in extracting the relevant, growing length scale 
from the simulated spin configurations. Here we calculate the two-point correlation function, G(r) — (s(ri)s(ri +r)), 
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FIG. 8: Simulation data for the correlation length £ (solid curves), together with three-parameter nonlinear fits (dashed 
curves), shown versus time t for lattices {4,4,60} (14,400 sites), {3,6,60} (10,800 sites), and {3,7,8} (11,173 sites). Linear 
scale (a) and log- log scale (b). In the legends, "raw" refers to data not corrected for thermal fluctuations, and "smoothed" 
refers to data with thermal fluctuations filtered out as discussed in the text. The filtered data give an estimate of the growth 
exponent of n « 0.3 for the Euclidean lattices, compatible with the theoretical value n = 1/3. However, the effective exponent 
obtained for the hyperbolic lattice is much lower: approximately 0.13. 

where r% is the position of lattice point i, and r, +r is the position of a lattice point a distance r away from i. Here, r 
is defined as the shortest path between two lattice points along the edges ("taxicab" or "Manhattan" distance). The 
correlation length £(£) is estimated as the first zero crossing of G(r) at time t. See Fig. [7j To reduce the effect on 
the estimate of thermal fluctuations in the bulk phases, we perform the simulations at relatively low temperatures, 
compared to T c . In calculating the correlation functions we also ignore isolated single spins and spin pairs |13j , which 
otherwise could distort G{r) for r = 1 and 2, as seen in Fig. [7] 



VI. NUMERICAL RESULTS 



The main numerical results of this study are summarized in Figs. [8] and [9j Figure [8] shows the length scale £ as a 
function of time for each of the three studied lattices, following a quench to T — 2.0. The data sets were averaged 
over 1,500 independent simulation runs for the Euclidean lattices, and 1000 runs for the hyperbolic lattice. Results 
are included, both based on the raw data, and with thermal fluctuations filtered out as discussed above. Significant 



differences are only observed for {4,4}, for which the quench temperature is not very far below T c . The growth 
exponents are here estimated by three-parameter, nonlinear fits to the time series, with the results given in the figure 
legends. For both the Euclidean lattices, the growth exponent comes out as n sw 0.3, consistent with the expected 
value of 1/3. For the hyperbolic {3, 7} lattice, however, the effective exponent is significantly lower, only about 0.13. 

In Fig. [9] we show results obtained by a different way of estimating the exponents. The data were divided into 
bins, each containing twice as many data points as the previous one ("octave binning"). We then performed a linear 
least-squares fit to log 10 £ versus log 10 1 in each successive pair of bins. Subtraction of a constant background was 
adjusted so that the estimated exponents became roughly independent of t. The resulting exponent estimates are 
seen to be consistent with those obtained by nonlinear fitting over the entire time interval. 

It is reasonable to ask whether the much lower effective growth exponent obtained for the hyperbolic lattice might 
be a result of finite-size saturation of the length scale. To check this possibility, we also performed simulations for 
smaller {3, 7} lattices with R between 3 and 6. As seen in Fig. 
10 6 MCSS, even for a system as small as R = 6. 
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FIG. 9: Estimates of the growth exponent n, based on the moving least-squares averaging method described in the text. The 
average values obtained are 0.30 for the two Euclidean lattices, and 0.12 for the hyperbolic {3,7} lattice. These values are 
consistent with the ones obtained by nonlinear fitting and shown in Fig. IS] dy refers to the constant background subtraction 
applied to the data to obtain exponent estimates approximately independent of t. 
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FIG. 10: Correlation length £ versus time for {3,7} lattices with R = 3,..., 6. These results indicate that the low value of 
the effective growth exponent observed with R = 8 (see Figs. [8] and |)J is genuine and not due to finite-size saturation of the 



domain size within the simulation time of 10 



8 
MCSS. 



VII. DISCUSSION 

Here we have reported a preliminary, numerical investigation of the dynamics of phase separation in a model with 
a single, conserved, scalar order parameter (Model B) confined to the vertices of a lattice obtained as a regular tiling 
of a hyperbolic plane. The results are compared with those obtained for the same model on two different Euclidean 
lattices. The latter show power-law domain growth with an exponent of approximately 0.3, near the theoretical result 
of 1/3, which is valid for the limiting case of an isotropic continuum system at zero temperature. Considering that 
the simulations were performed for anisotropic, discrete systems at a finite temperature, this agreement is convincing. 

However, for the hyperbolic {3, 7} lattice with R up to 8, we observe much slower growth, consistent with a power 
law with an effective exponent of about 0.13. Whether or not this is indeed power-law growth or something else, we 
leave open for future, theoretical investigation. The effect is possibly related to the fact that for large R, most of the 
spins are located near the free surface. It may also be related to the mean-field nature of the phase transition at T c 
in the hyperbolic case, which indicates the existence of effective long-range interactions. 

For the future we also leave a numerical investigation of phase ordering with non-conserved order parameter (Model 
A [9]) on hyperbolic lattices. In this case, the growth exponent in the Euclidean case is known to be 1/2. 
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